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C/3 ' In this paper we study the non-Gaussian features of the primordial fluctuations in loop quan- 



turn cosmology with the inverse volume corrections. The detailed analysis is performed in the 
single field slow-roll inflationary models. However, our results reflect the universal characteristics 
of bispectrum in loop quantum cosmology. The main corrections to the scalar bispectrum come 
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from two aspects: one is the modifications to the standard Bunch-Davies vacuum, the other is the 
corrections to the background dependent variables, such as slow-roll parameters. Our calculations 
show that the loop quantum corrections make /nl of the inflationary models increase 0.1%. More- 
over, we find that two new shapes of non-Gaussian signal arise, which we name T\ and J- 2- The 
former gives a unique loop quantum feature which is less correlated with the local, equilateral and 
single types, while the latter is highly correlated with the local one. 
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I. INTRODUCTION 



As a non-perturbative and background-independent theory, Loop Quantum Gravity 
(LQG) pj-yj has achieved great successes in past years: derivations of the quantized area and 
volume operators 4-7], calculations of black holes entropy jg] and Loop Quantum Cosmol- 



ogy(LQC) 9J, etc. And the nonperturbative quantization procedure of 



iQG is also valid for 



a more general class of four- dimensional metric theories of gravity [lGj-[l2| . As an example 
of LQG, LQC gives a quantization scheme of LQG for a symmetry-reduced model in the 
homogeneous and isotropic Friedmann-Lemaitre-Robertson- Walker universe. The discrete 
spacetime geometry in LQC scenarios predicts a non-singular bouncing universe in some 
simplified models, which satisfies most of the astronomical and cosmological observational 
constraints. Although the quantum correction effects are being diluted with the expansion 
of our universe, it remains present in a weaker form, especially on/near the super- horizon 
scales. 

Recently, the gauge invariant cosmological perturbation theory has been systematically 



constructed in 



13Ml5|| for inverse volume corrections and in 



18 



id 



19 



21 



171 for holonomy correc- 
22j . The inverse volume 



tions. Some relevant applications have been considered in 
corrections and the holonomy corrections are two main quantum corrections in LQC. The 
inverse volume corrections come from the quantization of the inverse of the volume operator 
in LQG. The inverse volumes exist in the Hamiltonian constraint of gravity and the usual 
matter Hamiltonian, especially in the kinematic terms. Since the volume can be taken the 
value zero, there does not exist well defined inverses of the volume operator. Fortunately, 



with the Thiemann trick 



23] , we can construct well defined inverse volume operators, which 



bring the quantum corrections. While the holonomy corrections arise from the loop quanti- 
zation based on the holonomies instead of the direct connection. The holonomy corrections 
become important when the energy scale of our Universe approaches the Planck one. Both 
the modifications to the scalar and tensor primordial power spectra from the inverse volume 
corrections are carefully investigated by the authors of 181 ] . Their results show that the in- 
verse volume corrections could give rise to the enhancement of the power spectra on the large 
scales, i.e., a red-tilt one. However, some other mechanisms such as the non- commutative 



geometry or string theory 



24 



26 



271 ]. could also lead to similar features. Therefore to seek 



for signature of loop quantum cosmology, the study for the non-Gaussianity features in loop 



2 



quantum cosmology is necessary. 

Because the primordial non-Gaussianities are quite helpful to distinguish inflationary 
models, so far a lot of papers have been devoted to studying the non-Gaussianities in different 



28 



291 ] . Inspired by the 



inflation models, see the relevant references in the nice reviews 
studies of 18j, in this paper we mainly consider the non-Gaussianities from the inverse 



volume corrections in LQC. The reason for consider inverse volume corrections only is as 
following. We denote <5; nv as the correction term coming from the inverse volume operator 
and <5hoi as the correction term from the holonomy corrections. We can estimate the inverse 
volume correction as 21 ] 

where the Planck density pp\ is assumed as the quantum gravity scale. From the above 
expression, we can see that the inverse volume corrections behave very differently from what 
is normally expected for quantum gravity. For low densities, the holonomy corrections is 
small, but the inverse volume one may still be large because they are magnified by the inverse 
of <5hoi- For an example, the small holonomy corrections of size 8^ \ < 10~ 6 then requires 
the inverse volume corrections larger than 5 inv > 10~ 6 even at scale p ~ 10~ 9 ppi. These 
novel features make the investigations on the inverse volume corrections more interesting 
than the former at sub-Planckian inflationary scales. So we only consider the inverse volume 
correction in this work. 

Explicitly, in the perturbation theory in LQC, the inverse volume operator can be cap- 
tured by a correction function such as a ~ 1 + « 5mv — 1 + «o( a mv/o) CT , where a is the scale 
factor of the FLRW universe and Oi nv is introduced to describe the characteristic scale of the 
inverse volume correction, which is not the Planck one in general. When a- imr /a < 1, we can 
ignore the correction term. However, if ai nv /a < 1 during inflation, one cannot neglect in- 
verse volume corrections. In this case, the correction approximates ao5- mv (k) ~ 5(ko)(ko/k) a , 
where k and fco are, respectively, the considered perturbation wave number and some char- 
acteristic number involved in the inverse volume correction. In addition, many works about 
LQC imply that a G [0, 6] From the form of the inverse volume correction, we can see 
that a small sigma corresponds to a small the inverse volume correct, vice versa. Therefore, 
as an example, following [ljj, we take a = 2 in this paper. For other a values the behavior 
will be similar. Furthermore, in terms of spherical multiples the wave number could be 
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expressed as k ps 10 _4 /i/, with h representing for the reduced Hubble parameter h = 0.7 and 
I for the spherical multiples. In the typical linear regime of Cosmic Microwave Background 
(CMB), the multiples I range in 2 < I < 1000 or more. Given all the mentioned observations, 
we expect some new features in the non-Gaussianities will arise. The purpose of this work 
is to investigate the characterized sizes and shapes of bispectra in LQC scenarios. 

Note that the bis pec tra for the single field slow-roll inflationary model have been calcu- 



lated in the papers 



30 



M34I]. For simplicity and comparison, we study the simplest single 
field slow-roll inflationary model in LQC. Our results show that the quantum corrections 
mainly come from the third order interaction Hamiltonian, the corrected vacuum state and 
the corrections to the slow-roll parameters. Our paper is organized as follows. In Sec. [TTJ 
the canonical formulism and the slow-roll inflationary model in LQC scenarios are briefly 
reviewed. In Sec. IIIH we study the power spectrum in LQC and recover the previous results. 
The effect of the inverse volume corrections on the non-Gaussianity in LQC is investigated 
in Sec. IIV1 The detailed analysis for the sizes, shapes and shape correlations is presented 
in Sec. HVknd Sec. [V] respectively. Throughout this paper we set 87r7G = 1 and Einstein's 
summing convention is always adopted. 



II. REVIEW OF LOOP QUANTUM COSMOLOGY 

The framework of LQG/LQC will be briefly presented in this section. Firstly, we discuss 
the canonical formalism in LQG, and then introduce the dynamics of slow-roll inflationary 
models in LQC scenarios. 



A. The canonical formalism in loop quantum gravity 

In the framework of LQG , the spatial metric as a canonical field is replaced by the 
densitized triad Ef, defined as 

E?:=\det(ei)\el (2) 

where e" is the inverse of the cotriad e\ related to the spatial metric by q a b = e^. The 
canonically conjugate variable to the densitized triad is the Ashtekar-Barbero connection 
A l a := + 7-fQ, where K l a is the extrinsic curvature and 7 ~ 0.274 is the Barbero-Immirzi 



4 



parameter |3 71. |38|. The densitized triad Ef and the Ashtekar-Barbero connection A l a satisfy 
the following commutator relation 

{Ai(x),E*(y)} = 6 b a 6i6 3 (x,y). (3) 

The spin connection T l a is defined such that it leaves the triad covariantly constant and has 



the explicit form 



r * = -e^Ad^eUik + ^e c k e l a d [c e%5 lm ) 



cj 



(4) 



In the new Ashtekar variables, the Einstein-Hilbert action can be expressed in the canonical 
form 



S 



EH 



dt 



where the Diffeomorphism constraint is 



%r,AN a ] = ~ / d 3 xN a (d a Ai - dtA^E'j - A{d h E) 
And correspondingly the Hamiltonian constraint can be expressed as 



The Gaussian constraint is 
<3 \A 



rpc rpd 

d 3 xNd h j k 



y/\detE\ 



2d e r i d + 4 m {r?r2-K?K2) 



d?xK\d a E« + e^TiEt + e%KiEt) 



(5) 



(6) 



(7) 



(8) 



which can be solved through standard procedure [l^j ]. Thus, solutions for the scalar mode 
perturbations are completely determined by the Hamiltonian constraint and the Diffeomor- 
phism constraint. 



B. Slow-roll inflationary models 

In this subsection, we shortly review the inflationary dynamics of the Friedmann- 
Lemaitre-Robertson- Walker Universe in LQC scenarios. The modified Friedmann equation, 
Raychaudhuri equation and Klein-Gordon equation are respectively [ljj] 



y?" + 2U 2 <p' (l-^f)+ vpV,M = 0. (H) 



2 



where T-L = f^, = du/dp, p = a 2 and a prime represents the derivative with respect to 
the conformal time, a and v are the correction functions for the inverse volumes and they 
read a ~ 1 + ctn&nv, P ~ 1 + ^ o^ii 



Following 19 



2 If . the slow- roll parameter e can be straightforwardly calculated as 



1 n ' 

= e (l + 7 Aw) , (12) 

where eo denotes for the usual slow-roll parameter. And the explicit form of coupling constant 
7 e reads 

le = So + a Q + v ( ^ - 1 J • (13) 



19]. So is determined 



2e 

Typically we set «o = 0.06, uq = 0.17 and eo = 0.01 in this paper 
by quantum correction and the analysis after Eq.( H3l) implies that it takes O(10~ 3 ). As 
mentioned in Introduction we set a G [0, 6] in this work. From above expression we can 
easily estimate that 7 e is of the order (9(10 -3 ). Formally we can also express the another 
slow-roll parameter j] as 

71 = 1 



Hip' ' 

= 77o(l+7rAnv) , (14) 

where 770 denotes for the usual slow-roll parameter as eo- 

As in the usual situation of the single field inflation model, we can assume that the two 
slow-roll parameters, e and rj , take roughly the same order as 10 -2 . Thus the typical values 
of coupling constant 7^ and 7 e are of the order C(10~ 3 ). The terms proportional to the Si QV 
represent for the inverse volume corrections. Here, we should emphasize that the subscript 
"inv" is introduced to avoid confusion with perturbations, such as S(p. 

III. POWER SPECTRUM 

In this section, we will firstly review the formalism of scalar perturbations in LQC; then, 
derive the second order Hamiltonian; and finally calculate the primordial power spectrum 
in the spatially flat gauge. 
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A. Formalism on the scalar modes 



Consider the scalar perturbations only, the general form of a perturbed metric around 
the isotropic FRW background is 

ds 2 = a 2 (r){ - (1 + 2<p)dr 2 + 2d a Bdrdx a + (1 - 2i))5 ab + 2d a d b E dx a dx b } , (15) 

where the scalar factor a is a function of the conformal time r, and (<f), ip, E, B) are the four 
scalar metric perturbations. In the perturbation theory, the triad can be described by 

E- = E? + 5E? , (16) 

where 

Ef = pSf, 5E? = -2p#? + p(6?A - d a d t )E. (17) 

The perturbed triad is described by the spatial part of the perturbed metric ip and E. Here 
A is the laplace operator in the flat space. Similarly, the perturbed lapse function and shift 
vector can be described by the other two scalar metric perturbation <fi and B respectively, 

5N = Ncj), 5N a = d a B. (18) 
The extrinsic curvature can be perturbed as 

K = Kl + 5KI = M a + SKI (19) 
For a general triad (fT6l) . the linearized spin connection becomes 

sr a = Y/« d > 5Eh r ( 2 °) 

As described above, the symplectic structure also splits into two parts, one for the back- 
ground variables and the other for the perturbations, 

{hp} = 3^, {5Ki(x),5E){y)} = , (21) 

where the background variables are defined by 

p= wJ EXrf3x ' k = w I K6td3x • (22) 

Here Vq is some artificial finite volume. 



In this paper, the matter part is represented by a scalar field <p. Similarly, we split the 
field tp and its conjugate momentum tt into homogeneous part and inhomogeneous one as 
well 



<p = p + Sip, 7i = 7i + 6n . (23) 

Hence, the basic Poisson brackets are reduced into 

{^,7f} = 1, {6<p(x),6Tr{y)} = S 3 (x-y). (24) 

For simplicity, we introduce the LQC formalism with finite cell Vq rather than the whole 
R 3 region in above description. But the unphysical feature of Vq can be remedied by lattice 
refinement model [3] . Since our followin g ca lculation only involves <5i nv we adopt the lattice 



refinement parametrization procedure in 



B. The second order Hamiltonian 



191 ] to eliminate the effect of artificial volume Vq. 



According to [14j, the quantum corrected second order Hamiltonian constraint can con- 
veniently be written as 

+ f d 3 x5N[ufiP+f)W] , (25) 
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where 



•'grav 



-4(1 + f)k^p5 c 3 5Ki - (1 + g)-=6i6E? + -=d c &bE] 

VP VP 



(26) 



(2) 
grav 



2k 





= (i + A) 




= o, 


4 1} 


= f /2 ((i 



^EI6EtS% + - (1 + h)^(dJ>E®(dJiEi) 

7l57T 

*2J 



p 



3/2 



2p 3 / 2 2p 



6>5Ej 



(2) 



-'7T 



Sit 2 

^ 1 + 9i)^-(l + 92) 



yr/ 2p 
TtSirSiSE? 1 vr 2 



(27) 

(28) 
(29) 
(30) 



p3/2 2 p 2p 3 / 2 



;i + 9s) 



8p 2 4p 2 



^(l + g 5 )Vp8 ab d a 8<pd b 8<p, 



(31) 
(32) 



(2) 



^3/2 



1 5 J '5i? c 
(1 + ge)-V >w (Cf>W + V„(<p)6<p-^ + 



(SlSEtf 5%5E^E( 



ip 2 



4p 



1 1- 



15] or in the Appendix B of 



(33) 



where the definitions of the counterterms can be found in 
And the perturbed second order diffeomorphism constraint can be expressed as 

d s x6N a [pd a (S d k 8K k ) - p(d k 5K k ) - k5 k (d d 5E d k ) + (ird a 6<pj\ • (34) 

Based on this corrected Hamiltonian, we can get the homogeneous and inhomogeneous part 
of matter field as 



vr = Cp'p/v, 5k = p{ V - ^(1 + h)<p] (1 - 9i) + 



(35) 



where fi and g\ are the counterterms. 



C. Power spectrum 

Calculations can be simplified greatly in the spatially flat gauge (if) = 0, E = 0), because 
the perturbed triad vanishes (SEf = 0) in this gauge. Here we fix the gauge after having 
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put quantum corrections in Hamiltonian and having checked consistency. In contrast, in 



references |39j, the authors fixed the gauge beforehand. We believe our treatment is more 
consistent. By solving the constraint equations, the perturbed lapse function and shift vector 
read 



AB = - 



(36) 



2VUl+f 
And the extrinsic curvature is 



' 3 {0V - v\i + h)<t> + fWAWv} ~ 3 ^(! + f)<t> ■ (37) 



a6K i a = -8\H{l + f)<j>-d a &B, 



(38) 



where ak = H. 

In the spatially flat gauge, the total second order Hamiltonian becomes 

~3pa a 2 pip' 4 



+ 
+ 



d 3 x 

s 

2 p2 ■'mr, 



2 , 2 ^^(1-^ + 2^-2/) 

v^M + *££vM(i + f 3 - f) s<p 



JL(1 - gi ) 5ip > 2 + p l(l + g,)5 ah d a 5 V d b 5 V ) 



(39) 



where 9 is also a correction function for the inverse volume and a 2 = 99. In this gauge, the 
dynamical inflaton perturbation 5<p coincides with the Sasaki-Mukhanov variable u = z(, 
with 



n L 



u )S 



(40) 



Then, one can derive the Mukhanov equation 19J 



u 



{c 2 s A + —)u = 

z 



(41) 



where c s is the propagation speed of the perturbation. The solution of the above equation 

is m 



u(k, t) 



H 



H 

7W 



-ikr 



-ikr 



1 + ikr 



X 



2 (7 + 1 



-{l + ikr)S h 



ikF [^) T + {k 2 T 2 ' 



(42) 
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where 



A 



-6 h 



C 



(43) 



^pjy and 5 = 5(k )/a . The latter 



2(<7+l, 

where \ — Wq(1 + a/6)/3 + «o(5 — <r/3)/2, C - 
variable 5(&o) i s constrained by the cosmic observational data [2l|, such as Cosmic Microwave 
Background and Large Scale Structures. For the specific inflationary models with a quadratic 
potential and a = 2, 5(k ) ~ (9(10~ 5 ). In this paper we take «o ~ (9(10 -2 ), i.e., the variable 
So is of the order O(10~ 3 ). Moreover, Eq. (|42p tells us that, on the one hand the inverse 
volume corrections become important for long wave modes with k <C k^ on the other 
hand long wave modes cross horizon earlier than the short ones. It means that the inverse 
volume correction will leave more hints on large scales than small ones. These features are 
much different from those of the inflationary models with higher derivative terms such as 
K-inflation in Einstein gravity. 

Using the canonical quantization, we have 



C(k, r) = C + + C = C(K r)a % + C*fe r)a f 



(44) 



where ((k, r) = u(k, r)/z(k). Then the two-point correlation functions of curvature pertur- 
bations can be calculated straight forwardly as 



\u\ 



(27T) 3 6(k + k'] 



H 2 



2e k 3 



{1 + 7 s £ inv } 



with 



7s = v (l + a/6) + aa /2e - xlip + 1) 



(45) 



(46) 



where the first and second terms in (1461) comes from the z(k) factor in the gauge transfor- 
mations, and the third term attributes to the modifications of vacuum state ( 1421) . 
Finally, we can get the primordial power spectrum of curvature perturbations as 



fc 3 



if 2 



2n 2 



This result is in agreement with that in 



47r 2 e 



fx 



21] . When all the corrections vanis 



(47) 
lis result 



is back to the case of the single field inflationary model in Einstein gravity |30H34| . The 
primordial power spectrum of curvature perturbations and angular power spectrum are 
plotted in Fig. (U where the dotted (purple), dashed (deep blue) and the solid (light blue) 
curves correspond to different values of the parameter C (C = 0,4 x 10~ 4 ,3 x 10~ 3 ) which 
was defined after Eq.(]4"3"l). 
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IV. BISPECTRUM 

In this section, we firstly derive the third order Hamiltonian in the spatially flat gauge; 
then calculate the three-point functions of the primordial curvature perturbations; and fi- 
nally figure out the sizes and shapes of the bispectrum. 



A. The third order Hamiltonian and In-In formalism 



We can get the corrected third order Hamiltonian 



»(3) 



(fx 



(48) 



where the expressions for fj^-L, i3y , are complicated, and we list them in Ap- 

pendix |A] 



12 



The perturbed third order diffeomorphism constraint is 



@®[8N a ] := - d 3 x5N a 

7 



2p 



\d a d c 5Et)5E b 3 + -fdJKiSE) - —^{dtdMDSE) 



^d b 8K{8E b - -ei k d c 8E c k d b 8E b - ^K{d b 8E) + 5nd a 5<p 
Zp 



35] 



We calculate the non-Gaussianity in the interaction picture 

(CV)> = (UrlC 3 (r)U int (r,To)), U int = e -'£ 
Up to the first order, we have 

(C 3 (r)) = -. f "rfr'([C 3 (r),^ nt (r')]}. 



(49) 



(50) 



(51) 



B. Sizes and shapes 



Based on the second order anomaly free perturbative LQC theory, we combine Eq.(j48J) 
with Eg. (1491) and arrive at the third order interaction Hamiltonian with some counterterms 



= ^( 3 ) + 0( 3 ) 



d 3 x 



^(1 - gi ) + y (1 + g 5 )5 ab dJ V d b 5^ 



+?(1 - g x )(PB)5<ffd a 8<p - ^(1 + 2/)K 2 3 
i/ a 



V 



+^-(l + 2A-^V + 



+— V tVW ,(<p)6(p 3 



11 + /l - -/ r / ,2 
= 0</9 



(52) 



where the first three terms are in the leading order under the slow-roll approximation. In 
the following calculations, we only consider these terms. Comparing Eq. ( |52l) with those 



in 



30 



-|34I|. we could attribute two kinds of modifications in the interaction Hamiltonian to 
the inverse volume corrections. One comes from the modification of vacuum state (l4"2l . the 
other from the background dependent coefficients, such as (p,6,a, g±, - • •). However, the 
non-Gaussian signatures from the latter are contaminated greatly by the cosmic variance. 
Hence, we ignore the modifications in v etc., when we calculate the parts of three-point 
functions directly from the In-In formalism. In another words, we ignore the quantum 
anomaly behavior for the Hamiltonian in this work. In order to demonstrate the reasons 



13 



more manifestly, we take the first term in ( 152|) as an example. 



drx 



2u 



gi)5(p' 2 (f) 



ax d(p dip 



2u - v ( - 

2 V3 



1 



5 inv 5(p' 2 5ip 



(53) 



where the first term in the second line appears in the usual form, while the second one 
contains <5 inv = S max r a . Note <5 inv depends on k, there is a corresponding k related to S max . 
When we put these terms into the formalism fl5Tj) and perform the time integral, we will 
obtain such corrections 

P 



K 3 



G(h, k 2 , h) + 



K = h + k 2 + h 



(54) 



where the dots denotes higher order corrections. From the above expression, we can see 
that such term peaks at k 3> K. Since larger k gives smaller 5pi, and quadruple is the 
lowest detectable mode in CMB, k corresponds to I = 2. So Eq. (I54jl peaks on the very 
low I ^ 2 region where the cosmic variance dominates over the signals (See Figure 1 and 



Figure 2 in 



2l||). Although the non-Gaussian features are presented on both large scale and 



small scale, above analysis implies that we can ignore the effects on small scales. Hence, 
we can ignore such terms in our following calculations. Once we ignore the corrections in 
the background dependent coefficients in Eq. (15211 . the form of the third order interaction 



Hamiltonian reduces to the usual one 



30 



34]. 



In order to eliminate the terms proportional to the linear equations, we need to do the 
field redefinition 



C = Cc 



1 - 

2 , n a-2//- n2 



(55) 



Then, the interaction Hamiltonian can be reduced into the simple form 



^int — 



According to 30h34|. after a field redefinition of the schematic form ( 
the correlation function will contain two terms 



(56) 



C c + X(c then 



(COsiXfoXfo)} = (Cc(xi)Cc(x 2 )C c (x 3 )) + 2A (C(x 1 )C(x 2 ))(C(x 1 )C(x 3 )) + cyclic , (57) 
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where the first term can be computed by the In-In formalism and the second term comes 
from the field redefinition ( = ( c + 
Let us firstly calculate the first term 



(Cc(t, fci)C c (r, k 2 )Cc(r, k 3 )) 



dr(0\ Cc(r,fci)Cc(r,A; 2 )C c (r,A;3),wS(r) |0) 



(27T) 3 S(£ k t )H* y? 4 4 E^-^ fc 

[nL* 2 (W]# 4 # 2 * 




(58) 



where z(k) and C are defined in Eqs. fl4"U|) and (143 p . respectively. Here, we emphasize again 
that corrections from the background dependent coefficients are neglected in the above 
expressions. However, in the next calculations of the parts from field redefinition, such 
corrections should be taken into account for consistence. Because the background dependent 
coefficients (</? /T-L\* , • • ■ ) take the value at the moment when the corresponding mode crosses 
horizon (r* ~ k' 1 ), they contain corrections such as (ko/ki) a . Of course, these terms also 
peak at the points where k ^> ki, they become important, particularly, in the squeezed 
triangle limit {k\ <C k 2) k 3 ). 

The contributions from field redefinition can be decomposed into two parts, one is 



(2vr) 3 5(^£ t ) 



AH 4 



[nti^w] 



i=l 



(59) 



and the other is 



2H A 
[UU^k^kf] 



Y,C 2 {k t )z 2 {k l )k l k 2 j 

Art 



(60) 



where H = ^ and the overdot stands for derivative with respect to the cosmic time. The 
Ci(ki), C 2 (kj) terms can be read from Eq. fl55|) by substituting r with fcr . 

In summary, we conclude that the forms of interaction Hamiltonian are exactly the same 



as the usual one 



30 



34l | . however, the inverse volume corrections 5 mv will make some con- 



tribution to the bispectrum. There are mainly two sources, one is the modifications to the 
standard Bunch-Davies vacuum, the other comes from the z(k) factor in the gauge trans- 
formation ((k) = u(k)/z(k). Furthermore, we can expand the above results in terms of 
<W = 8o(ko/kiY an d obtain 



(27r) 3 5(^ ki) jingle {h, fc 2) h) + k 2 , fa) + Fi(fa, fa, k 3 ) 



(61) 
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where 



single 



nL(2* 



( 3eo - 2 V0 ) kf + e J2 hk) + 86 Ei>i ^ 

i i^j 



K 



(62) 



is the usual leading term and 

3 



•F 2 



Z=l 



ki 



single 



(63) 



(2tt) 4 P c 2 



nL(2^; 



1 ' #3 1 



(64) 



are the inverse volume correction terms, and here the relevant coefficients are 



2 

C z 

c 



2C-C Z , 



2C . 
2C 



-S 



2e7 



2(a + l) 



(65) 
(66) 
(67) 
(68) 

(69) 

terms in the denomi- 



Particularly, we figure out that the corrections from 
nators in Eqs.( l58p . ( 15 9 p and (160 p are absorbed into T\ shape, while all other corrections are 
collected in T%. In the next section, we will find T\ shape provides an unique signal from 
LQC mechanism, and more importantly, this signal is independent of the inflationary mod- 



els, because 



[\^ =1 z 2 (hi) terms always appear in the gauge transformations ( 140]) . Namely 



it is an universal signal in the LQC scenario. 

The single shape is the usual one, while T\ and arise only in LQC scenarios. Sizes 
of the two new parts are proportional to parameters {u\ ,•■•), which are in the order of 
O(10~ 3 , 10~ 4 ). That is to say that, these new non-Gaussian features from LQC are smaller 
than those of usual inflationary models in Einstein gravity by a factor at least 10~ 3 . Although 
this does be a tiny number, considering these features sourced by quantum effect, this factor 
is not small as initially expected. Furthermore, as stated above, because the inverse volume 
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FIG. 2: Single (a) and F\ (b) shapes. The z-axis is x\x\T{x%^ Xs)/P^, here we set x\ = 1. The 
slow-roll parameters take the typical values eo = % = 0.01. Here and hereafter we set a = 2. 

corrections in the interaction Hamiltonian can be neglected, i.e. we could use the usual 
Hamiltonian directly in the In- In formalism by substituting the Bunch-Davies vacuum state 
with the one given in Eq.( )42]) . we argue that the expectations on the sizes of bispectrum 
should hold for any inflationary models in LQC scenario. 

In Fig. H] (a) (b) and Fig. [3] (c), we plot x\x\x\T{x\, x 2 , x 3 )/P^, with Xi = ki/k\. 
The difference between shapes T\ and Ti is plotted in Fig. [3] (d). We can see that all 
the three shapes peak at the squeezed limit (x 2 = l,x 3 = 0), however, the substructures 
are different among them. Compared to the single shape, T\ shape upraises at the corner 
(x2 = 0,X3 = 0.5), while T% flattens at the same point. From Fig. [3] (d), we can also see 
that J-2 peaks more dramatically in the squeezed corner than T\. 



V. SHAPE CORRELATIONS 

In the Figures [2] and [31 all the three shapes looks similar. In order to figure out the 
differences among shapes more quantitatively, we need to calculate the correlations between 
them. 
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10 X 3 i.o X 3 

FIG. 3: T 2 (c) shape and T2 — J 7 ! (d). The coefficients ,uj' 2 ,uj 3 ) in T\,Ti are set unity. 



Firstly, we define a 3D shape function 40h42| 



S(h, k 2 , k 3 ) = ^{k.hhfHh, h, h) , (70) 

where iV is a normalization factor which will not affect the following calculations. 
Then, we construct the products of two shape functions 

F(S, S') = [ S(ki, k 2 , h)S'(h, k 2 , ksMh, k 2 , k 3 )dV k , (71) 

where u(ki, k 2 , k 3 ) is a weight function and Vk is the integration domain constrained by the 
triangle inequality. 

Finally, we arrive at the 3D shape function correlator 

£(g g'\ _ S ) 

y/F(S, S)F(S', S>) ' 

This quantity describe the cross correlations between two different shapes. 

Furthermore, for a large series of well-motivated shapes, the above descriptions can 
be simplified. One can define the distance described by k from the origin of (^1,^2,^3)- 
momentum space to the particular triangle slice which is perpendicular to (1, 1, 1) direction, 

k = -{hi + k 2 + k 3 ) , (73) 
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then, we introduce another two new variables 



h = k(l-P), (74) 
k 2 = h{l + a + (3) , (75) 

h = \k{l -a + f3). (76) 

In the domain constrained by the triangle inequality, 0</c<oo,0</3<1 and —(1 — /3) < 
a < 1 — (3. For the classes of models with homogeneous shape, which means that the powers 
of wavenumber in shapes are homogeneous, the k dependence in the 3D shape function can 
be separated 

S(k u k 2: k 3 ) = f{k)y{a, (3) , dV k = dk x dk 2 dk z = k 2 dkdad(3 . (77) 
In fact, for the models considered in 
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42j] and here, f(k) oc const. Hence, we can focus 
on the 2D shape function S^(a,(3) and the integral over k cancels when one calculates the 
shape function correlators. 
Further, one can introduce 

a = l-x, fi = xy/3, (0<x,y<l), (78) 

to square the integration regime. By using these variables, the integral measurement becomes 

dadf3 = xdxdy . (79) 



From the above expression, we can read t 



re weight function w(x,y) = x. This choice works 



well for all the shapes mentioned in 4CH42j . however, for our new shapes T\ and T 2 the 
correlation matrices c mn defined in ( )84|) . suffers from divergence. In our calculation we 
therefore use a new variable £ 2 = x 

dad(3 = t, 3 d£dy , (80) 

to eliminate such divergence. Using variables (£, y), we can decompose the shape J^(£, y) on 
any triangle slice with analogous radial polynomials R m {0 and shifted Legendre polynomials 

P n (y) 

y(^y) = Y, c ^RUOPn(y) , (81) 
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where the first few Rm(Q eigenfunctions are 



Ro = V2 , R 1 = y/4(-2 + 30 , R 2 = V6(3 - 12f + 10£ 2 ) 
R 3 = v/8(-4 + 30£ - 60£ 2 + 35£ 3 ) 



(82) 



And P n (y) eigenfunctions are 



P = l, Pi = V3(-l + 2y) , P 2 = V5(l - 6j/ + 6j/ 2 ) 
P 3 = V7(-l + 12y - 30y 2 + 20y 3 ) 



(83) 



Thus, one can define the correlation matrix c r 



dC f d y e^{U)R m ^)Pn{y) ■ 

o Jo 



(84) 



The c mn matrices for local, equilateral, single field model, T\ and T 2 are listed in (185]) and 
Fig. HI The explicit definitions of 2D shape function S^{^,y) can be found in Appendix IB1 



/ 1.00 -0.14 0.03 0.00 \ 
0.38 -0.07 0.02 0.00 
0.04 -0.01 0.01 0.00 
\0.02 0.00 0.00 0.00 J 

( 1.00 -0.19 0.04 0.00 \ 

-0.26 0.09 -0.01 0.00 

0.44 -0.10 0.02 0.00 

\ -0.32 0.08 -0.01 0.00 / 



/ 1.00 0.40 -0.12 0.01 \ 

0.68 0.27 -0.09 0.01 

0.21 0.07 -0.03 0.00 

\0.01 0.00 -0.01 0.00 / 



/ 1.00 -0.07 0.01 0.00 \ 

0.43 -0.02 0.01 0.00 

0.07 0.00 0.00 0.00 

\0.02 0.00 0.00 0.00 / 

/ 1.00 -0.12 0.24 0.00 \ 

0.30 -0.03 0.01 0.00 

0.07 -0.01 0.01 0.00 

\0.02 0.00 0.00 0.00 / 



■5) 



From Fig. H] we can see that T\ term differs from others explicitly, while T 2 is almost 
indistinguishable with the local form visually. 

Armed with the above results, one can calculate 2D shape correlator 



where the product is defined through 



(86) 



J Si, ™ ~ 



(87) 
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1 2 3 4 




FIG. 4: Correlation matrices in (|85p . Light blue denotes for +1, and the deep for -1. 





Local 


Equi 


Single 


Ti 


T* 


Local 


1 


0.826 


0.996 


0.722 


0.996 


Equil 




1 


0.874 


0.535 


0.826 


Single 






1 


0.707 


0.992 


?i 








1 


0.769 


?2 










1 



TABLE I: 2D shape correlators. 



The numerical results of the 2D correlators are listed in Tablet HI from which we can find 
that the correlations of shape T\ are low with all other four shapes, while shape Ti possesses 
high correlations with local (particularly), equilateral as well as single form. That is to say 
that J-2 is almost indistinguishable with local form, while T\ does be the unique signal of 
LQC. As already argued in the previous section, T\ is an universal shape in LQC scenario, 
so we can identify this shape as a new window for LQC scenario. 

Finally, let us estimate the parameter /nl- Here we focus on the two new shapes J-1,2, from 
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Tab. [J we can see that T% is highly correlated with the local form, while T\ is less correlated 

•loca 
NL 



with them. The contributions to f^ al from these two shapes can be easily estimated as 



A^° L cal ~ u x ^(^ local ,^ li2 ) x /^(single) ~ O(10~ 3 ) /^(single) , 

where u represents in short parameters (ui,uj' 2 ,uj 3 ) in fl63|) and flMl) . whose typical values 
are of order 0(lCr 3 ). Because f^l 3,1 (single) in the usual case is of order ~ e |43j, the 
contributions from inverse volume corrections is completely negligible. However, as argued 
before, our results should be robust for other inflationary models in LQC scenarios, especially 
those with large non-Gaussianities, such as K-inflation, DBI inflation etc. So, one can 
anticipate that in those models with large non-Gaussianities the features from inverse volume 
corrections in LQC scenarios might be observed. 



VI. CONCLUSIONS 



In this paper we investigated the contributions to the cosmic primordial scalar bispec- 
trum from the inverse volume corrections in LQC scenarios. We derived the interaction 
Hamiltonian, however, we found that the new interactions contribute greatly to the modes 
with k\ + k 2 + k 3 <C k. Because the scales corresponding to k is very large, here we take 
k ~ 0.00014Mpc _1 corresponding to quadruple mode I = 2, it means that the three wave 
lengthes in the bispectrum are all on the super-horizon scales. On so large scales the cosmic 
variance usually dominates over the signals, we hence neglected these new interactions in 
our calculations. That is to say, the interaction Hamiltonian we used is the same form as 
the usual one for the single field inflation models in Einstein gravity. This greatly simplifies 
our calculations, and more importantly, makes our results robust, i.e., our results should 
hold for other inflationary models in LQC scenarios. 

Although the Hamiltonian shares the same forms as the one in Einstein gravity, the inverse 
volume corrections 5- mv oc (fco/fc) CT still contribute to the bispectrum. Roughly speaking, there 
are two aspects, one comes from the deviations from the standard Bunch-Davies vacuum; 
the other attributes to the non-trivial gauge transformations ((k) = u(k)/z(k) from the 
spatially flat gauge to the observable curvature perturbations. Consequently, we obtained 
the three-point functions of the gauge invariant curvature perturbations. We found that, 
except for the usual single component in slow-roll inflation models, two new shapes arise due 
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to the corrections, namely T\ % 2- Furthermore, we performed a careful analysis on the new 
shapes. We found that, the whole profiles for all the three shapes (single, are visually 
similar, i.e. they peak at the squeezed limit. However, the substructures among them are 
different. Compared to the single shape, T\ shape upraises at another corner (See. Fig. [2] 
and [3]), while Ti flattens at the same point, and Ti peaks more dramatically in the squeezed 
corner than T\. In addition, we investigated the correlations among five shapes, including 
local, equilateral, single and T\ t i- The results show that Ti is highly correlated with the 
local type, while T\ is less. It means that the latter can provide a new window for probing 
the loop quantum mechanisms using cosmic primordial bispectrum information. Finally, we 
estimated the order of observable parameter Af^ 8,1 ~ O(10~ 3 ) x /^l (inflation) from the 
inverse volume corrections. 

The non-Gaussianity from the inverse volume corrections in LQC scenarios is tiny and still 
undetectable currently, however, considering they are generated by the quantum effect, which 
is naively expected of the order (9(GUT/Planck) fT ~ (9(10 _5 ) <T , our finding becomes non- 
trivial. Especially, the results obtained in this work should be generalized directly to other 
inflationary models with large non-Gaussianities, in which the inverse volume corrections also 
becomes large therein. In addition, in this paper we only investigated the non-Gaussianities 
from the inverse volume corrections, while ignored those from the holonomy corrections, 
in which the sound speed of scalar perturbations are typically changed. Besides, we only 
investigated in this work the bispectrum from scalar modes, left those from tensor mode 
unexplored. These topics are worth investigating further. 
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Appendix A: Interacting Hamiltonian 



In this appendix, we derive the perturbed Hamiltonian density and the perturbed dif 



feomorphism constraint up to the third order. According to 
includes two parts 



15[, the classical Hamiltonian 



3tf[N] = J? glav [N] + Matter [N] 



(Al) 



The gravitational Hamiltonian can be expressed in terms of the extrinsic curvature 
1 



•Kr^iN] = - / <PxN$) = - f d 3 xNei k -!¥k- 
and the matter part of the Hamiltonian is 



+ , (A2) 



^matter [N] = j d*xNfa + ij V + 



where 



fin 



* 2 ft EfEtdgipdbip . 



The diffeomorphism constraint is 



%rav[N a ] := /J 3 



(d a Al - d b A^)E" - Aid h E\ 



and the matter contribution is 



(A3) 



(A4) 



(A5) 



\N a ] := j d 6 xN a -nd aV . 



(A6) 



In order to obtain the perturbed Hamiltonian density and the perturbed diffeomorphism 
constraint up to the third order, we need the following two relations. We expand (det E) 5 
and (det E)~* to the third order as follows. 



(det E) 2 = p 



2p 



^ k JEl){8l8i8Et8E) 



12p 



+ 



1 



12p 



-8l8l8E?8E»8E c k + ... 



(A7) 
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and 



(detE)- 



3 

P 2 



ip 6 



2p 8p 2 
-(5 k c 5EtMV a 5Et5E") 



a\2 



" -SisMSEfSE^Et 



Up 



3 u c u b u a L 



; 5l5l5 k JE?5E b JE? + 



12 p3 b c a 

Thus the third order gravitational Hamiltonian density can be written as 



(A8) 



(A9) 



where 



fli- 



cs) 



jk j \ 2 d c r d 



2e 



jk 



VdetE 
6Ef6E$ 

VdetE 



(3) 



djr d + 2e> k 5E^E d k [(detE)-^d c ST d 



+2ei fc -Ej^-Efc[(det £ , )~^ 1 ^(9 c 5r^ + 2^ k EjE k [(det E)~^]^d c 5T d 

%5 kl 5 d 5Ef5E c j (d c d e 5El) + Xs^SEtS^SE^d^SEl) 
p2 ipi 



4p 



2p; 



(A10) 



Jk j k j -pm-pr 
1 ran c a 



VdetE 



(3) 



VdetE 



E c 5E d 

rpmrpn i J fc j XT™ %T^n 



+^*^[(det£?)- 



A -5 

4^2 



2p 



-5l5 h i5EC(d b 5E b k )(d a 8E?) 



(All) 
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, h 8E c 4 8E d . _ ^E^8E d ■ 

-4 \/^ e ™ K ? 5K a - <^==^JK™8K n d 

-ef5E^5E d k [(det E)^f^ mn K^K n d - ef E°E d [(det ^f^SK^SIQ 
-4 k 6E*E d k [(det Ey^V mn 5K™K n d - ef 5£^[(det Ey^e^ mn K™8K n d 
-4 k E'8E d k [(det E)-^ l ^ mn K n d 5K™ - 4 k E*8E d [(det Ey^e^ mn K™8K n d 
-4 k E d 5E*[(det E)^fV mn K^K n d - 4 k E°6E d [(det E)^f^ mn K™K" d 
-4 k E]E d k [(det EY^^JK-K- - ef £j^[(det Ey^^K^SlQ 

-etEptKdetEr^e^K- 

2k 2k 2 

-^8 k 8E d k 8E^ c + -8 d 8E^8E d k 8K k - ^8 d n 8K n d 8E]8K{ 

p2 p2 p2 

+^ T 8 d m 8E< j 8K™8K d + ^{8{8E^ - ^8 d 8 k 8 l a ^8E d 8E^ 

P 2 2p 2 2p 2 

1 1 2^ 

+^(^iC) 2 ^£? - ^^i^iC^f + -8E c 3 8Kl8 l a 8E? 

2p2 2p2 p2 

+\{^f{8 d n 8K n d ) - Xvj d c 8E<8K n d 8 l a 8E? - \K(&E<f 

p2 p2 A p2 

k 2 k k 

--(8 l b S l a SE?8Et)(8i8EC) - -^(8 l a 8En 2 (8 d 8K n d ) - -{8 l b 8l8E?8E*){8 d n 8K d 

P 2 2p 2 p " 

k 2 ^k 2 



2 



2 



Bp 2 4p 

k2 -8l8 k 8i8E?8E»8E c k + ^8l8l8 k 8E?8E>8E c k . (A12) 



2p2 2p 
The matter hamiltonian can be expressed as 

= ^vr 2 [(det + 7r57r[(det £)^]( 2 ) + ivr 2 [(det £)~¥ 3) 



5 I ' _3 
4^2 p2 



| 8l8j8E*8E b 3 
8p 2 Ap 2 



7T 



2 



2p2 



(W) 3 (ggjg^gKgg) 8l8 k h 8l8E*8E)8Ej 8l8l8 k a 8E«8E)8Ej , 

48 p3 g p3 12 p3 12 p3 J' ^ ^ 
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(3) 5 ij EpE?d a Stpd b 5ip ^ 8 ij E? 8E b d a 8(pd b 8(p ^ 8 ij EfE] 



+ 



2 v / detE 2v / det £ 

8 ai 8E b i d a 8yd b 8ip 8 k 8E c k 8 ab d a 8<pd b 8<p 



+ 



[(det£?)-3]W9 tt ^^ 



P 5 



(A14) 



fl® = ^Vd^F^)V + ~[(det£)^^ 



+ 



(3) 

(detS)s 

1 



i X Z?a\2 



8j8l8E?8E b 

4p2 



V, v (<p)6<p + p* 



{8 k c 8Ejf _ {5 k c 5Et){8l8j8E*5Ef 
48^3 



8p 3 



V c b k h blbE1bE)bEi 5l5t5 k JE?5E b 3 5Et 



V{<p). 



(A15) 



12p 3 12p 3 
Combined with the results obtained by [l^ . we get the perturbed Hamiltonian as follows. 



where 



and 



(2) 

grav 



and the matter Hamiltonian reads 



Walter = / 



4 3) +4 3) +4 3) 



(A16) 



(A17) 



(A18) 



(A19) 



d 3 x8N 



(A20) 



Here we have ignored the high order correction terms caused by the inverse volume, such 
as a because in the in-in formulism, these terms do not contribute to the non- 

Gaussianity. 

The diffeomorphism constraint up to the third order is 



- / d 3 xSN a 
7 is 



^ b j {d a d c 5El)5E) + iMKiSEj - y a \d b d c 5Et)5E) 



ydbSKlStf - —^ a k d c 8E c k d b 8E b - ^SKi^SE. 

Zp 



(A21) 
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matter 



[5N a ] := / d d x8N a 57rd a 5tp 



(A22) 



Appendix B: Definitions of shape functions 



For all the shapes considered in this paper f(k) = const., so we have 



S{h,k 2 ,k 3 ) = y(t,y). 



For the local shape 



for the equilateral shape 



5, 



local 



iEr-3 _i_ 7^3 _i_ L.3 
ru-^ nj<2 \ rug 

3hk 2 k 3 



S, 



cquil 



(fci + k 2 - h){k 2 + k 3 - fci)(fc 3 + fci - fc 2 ) 



for the single shape 

'-'single 

for the J- 1 shape 



(36 - 2 m ) kf + 6q fc^ + 86 fcffi/if 



x 5. 



single i 



and for the T 2 shape 



S 



hk 2 k 3 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 
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